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I present a brief overview of a variety of computational tools for su- 
persymmetry calculations, including: spectrum generators, cross section 
and branching fraction calculators, low energy constraints, general pur- 
pose event generators, matrix element event generators, SUSY dark mat- 
ter codes, parameter extraction codes and Les Houches interface tools. 

1.1. Introduction 

The Standard Model (SM) of particle physics provides an excellent descrip- 
tion of almost all physical processes as measured in terrestrial experiments, 
and is rightly regarded as the crowning achievement of many decades of 
experimental and theoretical work in elementary particle physics. 1 

As exciting as this is, it is also apparent that the SM cannot account 
for a wide assortment of astrophysical data, including neutrino oscillations, 
the matter-anti-matter content of the universe, the presence of dark energy 
and the presence of dark matter in the universe, and it doesn't include 
gravitation. Even before these astrophysical anomalies became evident, 
it was apparent on theoretical grounds, mainly associated with quadratic 
divergences in the scalar (Higgs) sector, that the SM was to be regarded as 
an effective theory valid only at the energy scale of ~ 100 GeV and below. 
At higher energies, it seemed likely that some new physics must arise, which 
would be associated with the mechanism for electroweak symmetry breaking 

While a vast array of physics theories beyond the SM have been pro- 
posed, the general class of theories including weak scale supersymmetry 
seem to most naturally solve the theoretical ills of the SM, while at the 
same time they receive support from a variety of precision experimental 
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measurements. 2 Most impressive is the measured values of three SM gauge 
couplings at the weak scale: when extrapolated to high energies using the 
renormalization group group equations, the gauge couplings very nearly 
meet at a point under supersymmctric standard model evolution, while 
they miss badly under SM evolution. 3 Gauge coupling unification suggests 
that physics at scales Mqut ~2x 10 16 GeV is described by a supersym- 
metric grand unified theory, and that below Mqut, the correct effective 
field theory is the Minimal Supersymmetric Standard Model (MSSM), or 
the MSSM plus gauge singlets (since gauge singlets don't affect the running 
of gauge couplings at one loop). 

Supersymmetric models predict the existence of a whole new class of 
matter states at or around the weak scale: the so-called super-partners. 
Gluinos, charginos, neutralinos, squarks, sleptons plus additional Higgs 
scalars (H, A and H ± ) should all be present in addition to the usual states 
of matter present in the SM. In order to fully test the hypothesis of weak 
scale supersymmetry, it seems necessary to actually produce at least some 
of the superpartners at high energy collider experiments, and to measure 
many of their properties (mass, spin, coupling strengths and mixing), in 
order to verify that any new physics signal indeed comes from superpartner 
production. In addition, the properties of the superpartners will be key to 
understanding the next level of understanding in the laws of physics, per- 
haps opening windows to the physics of grand unification and even string 
theory. 

The key link between theoretical musings about various theories of 
SUSY or other new physics, and the experimental observation of particle 
tracks and calorimeter depositions in collider detectors is the event gener- 
ator program. Given some theory of new physics, which usually predicts 
the existence of new matter states or new interactions, the event generator 
program allows us to compute how such a theory would manifest itself at 
high energy colliding beam experiments. Thus, event generator programs 
function as a sort of beacon, showing the way to finding new physics in a 
vast assortment of collider data. 

Searches for new matter states at the CERN LEP2 collider have found 
no firm new physics signals. We thus conclude that the SM Higgs boson 
must have mass mH SM ~ 114 GeV, while the charginos of supersymmetry 
must have mass ~ 103.5 GeV. The Fermilab Tevatron is probing spar- 
ticle masses such as the gluino up to the 300 GeV level. The CERN LHC 
is a pp collider which is just now beginning to explore the energy regime 
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where the SM breaks down, and where new physics ought to lie. LHC is 
expected to operate at energy scales y/s = 7—14 TeV; this ought to be 
sufficient to either produce some superpartners, or rule out most particle 
physics models which include weak scale supersymmetry. 

As we enter the LHC era, it is important to review the available calcula- 
tional tools which are available, that aid in connecting theory to experiment. 
In this chapter, I present a brief overview of some of the publicly available 
tools. In Section [Oil I examine the various 

• sparticle mass spectrum calculators, 

and the Les Houche Accord files which provide a handy interface between 
these and event generator programs. Sec. 11.31 lists some 

• codes which calculate sparticle production rates, decay widths and 
branching fractions. 

Sec. 11.41 reviews 

• event generators for SUSY processes, including 

— multi-purpose generators, complete with parton showers, 
hadronization and underlying events, and 

— more specialized matrix element generators, which tend to 
focus on specific reactions. 

The Les Houche Event files allow parton level collider events to be easily 
passed into general purpose event generators so that showering, hadroniza- 
tion and underlying events can be included. Sec. II. 51 lists 

• codes relevant to supersymmetric dark matter calculations, 
while Sec. 11.61 examines 

• codes designed to extract fundamental theory parameters from sets 
of experimental measurements. 

The supersymmetry parameter analysis (SPA) project seeks to develop a 
uniform set of conventions which would allow unambiguous extraction of 
high energy model parameters from various collider measurements of su- 
persymmetric production and decay reactions. 

I note here that this Chapter is an updated version of the 1997 version 
by H. Baer and S. Mrenna which appeared in the volume Perspectives on 
Supersymmetry, edited by G. Kane (World Scientific). 4 
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1.2. SUSY spectrum calculators 

The first step in connecting supersymmetric theory to experiment is to 
begin with a supersymmetric model, and calculate the expected spectrum 
of superpartner and Higgs boson masses and couplings. The models we will 
be focusing on are 4-d supersymmetric quantum field theories with softly 
broken supersymmetry at the TeV scale. These models might be the low 
energy effective theory resulting from some even more encompassing theory, 
such as superstring theory, or a particular SUSY GUT model, or which may 
invoke some specific mechanism for supersymmetry breaking. 

The effective field theory is specified 5 by adopting 1. the gauge sym- 
metry, 2. the (super-) field content and 3. the Lagrangian. In the case 
of supersymmetric theories, the Lagrangian is derived from the more fun- 
damental superpotential and Kahlcr potential, and for non-renormalizablc 
models, the gauge kinetic function. The effects of SUSY breaking are en- 
coded in the Lagrangian soft SUSY breaking (SSB) terms. One must also 
specify the energy scale at which the effective theory and Lagrangian pa- 
rameters are valid. Since collider experiments will be testing physics at 
the weak scale Q ~ 1 TeV, while the Lagrangian parameters are frequently 
specified at much higher scales {e.g. Mqut or Mp), the renormalization 
group equations (RGEs) must be used to connect the disparate scales in the 
model. 

Once the Lagrangian parameters are known at the weak scale, then the 
physical (s)particle masses must be identified, often by diagonalizing the 
relevant moss matrices. Higher order perturbative corrections to the mass 
eigenstates- at least at 1-loop- are nowadays necessary to gain sufficient 
accuracy in the predictions. 6 

Numerous researchers have developed private codes to calculate sparti- 
cle masses given high scale model inputs. Here, we will focus only on pub- 
licly available codes, since these are available to the general user, and are 
frequently kept up-to-date and user friendly. The first of the publicly avail- 
able spectrum calculator codes to appear was the ISASUGRA 7 subprogram 
of the event generator Isajet, 8 in 1994. This was followed by SuSpect 9 
(1997), Softsusy 10 (2002) and Spheno 11 (2003). 

1.2.1. Isasusy, Isasugra and Isajet 

ISASUSY is a subprogram of the Isajet event generator which calculates 
sparticle mass spectra given a set of 24 SSB input parameters at the weak 
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scale. The program includes full 1-loop corrections to all sparticle masses. 
For Higgs masses and couplings, the full 1-loop effective potential is min- 
imized at an optimized scale choice which accounts for leading 2-loop ef- 
fects. 12 Yukawa couplings which are necessary for the loop calculations are 
evaluated using simple SM running mass expressions. 

The ISASUGRA 7 program starts with models defined at a much higher 
mass scale (e.g. Q = Mqut), and calculates the weak scale SUSY param- 
eters via the full set of 2-loop RGEs. 13 An iterative approach to solving 
the RGEs is employed, since weak scale threshold corrections which de- 
pend on the entire SUSY mass spectrum are included. Once convergence 
is achieved, then the complete set of 1-loop corrected sparticle and Higgs 
masses are computed as in ISASUSY . Since ISASUGRA employs full 2-loop 
running of gauge and Yukawa couplings including threshold corrections- 
while ISASUSY does not- the sparticle masses will differ between ISASUSY 
and ISASUGRA even for the same weak scale parameter inputs. 

A listing of pre-programmed ISASUGRA models include the following: 

• mSUGRA (or CMSSM) model: 4 parameters plus a sign 
(mo, mi/2, A , tan/3, sign(u)), 

• minimal gauge- mediated SUSY breaking (mGMSB, 4 param's plus 
sign plus C grav ) and non- minimal GMSB, 

• non-universal supergravity (19 param's plus sign) 

— SSB terms can be assigned at any intermediate scale M wea k < 
Q < M GUT , 

— non-universal Higgs model with weak scale n and vtia inputs 
in lieu of m Hu and m H , 

• mSUGRA or NUSUGRA plus right-hand neutrino (RHN), 

• minimal and non- minimal anomaly mediation (AMSB), 

• mixed moduli- AMSB (mirage mediation), 

• hypercharged AMSB. 

The related program RGEFLAV computes the complete flavor matrix 
structure of the SSB terms and Yukawa couplings, including CP-violating 
effects. 14 The webpage is located at http://www.nhn.ou.edu/~isajet/. 

1.2.2. Suspect 

SuSpect 9 runs the 2-loop MSSM RGEs to determine weak scale 
SUSY parameters in the mSUGRA, GMSB and AMSB models, 
and in the pMSSM (a more general MSSM model). One-loop 
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sparticle mass corrections are included. Some two loop correc- 
tions to Higgs masses are included. The webpage is located at 
http : //www . lpta . univ-montp2 . fr/users/kneur /Suspect/. 

1.2.3. SoftSUSY 

Softsusy 10 is a C + + code that calculates 2-loop MSSM RGEs to deter- 
mine weak scale SUSY parameters in the mSUGRA, mGMSB and mAMSB 
models, and in the general MSSM. i?-parity violating effects are possible. 
One-loop sparticle mass corrections are included. Some two loop corrections 
to Higgs masses are included. SoftSUSY calculates the complete flavor ma- 
trix structure of the MSSM soft terms and Yukawa couplings. The webpage 
is located at http://projects.hepforge.org/softsusy/. 

1.2.4. Spheno 

Spheno 11 is a Fortran 90 code that calculates 2-loop MSSM 
RGEs to determine weak scale SUSY parameters in the mSUGRA, 
mGMSB and mAMSB models, and in the general MSSM. One- 
loop sparticle mass corrections are included. Some two loop cor- 
rections to Higgs masses are included. The webpage is located at 
http : //www . physik . uni-wuerzburg . de/~porod/SPheno . html. 

1.2.5. Les Houches Accord (LHA) files 

A standard input/output file under the name of Les Houches Accord (LHA) 
files has been created. All of the above codes can create LHA output files. 
The advantage of LHA output files is that various event generator and dark 
matter codes (see below) can use these as inputs for generating collider 
events or dark matter observables. The specific form of the LHA files is 
presented in Ref. 15 

In addition, the ISASUGRA and Isasusy codes output to a special 
IsaWIG file, which is created expressly for input of sparticle mass spectra 
and decay branching fractions to the event generator Herwig . 

1.2.6. Comparison of spectra generator codes 

Several papers have been written comparing the SUSY spectra codes, 16 
although these tend to be all dated material, as the codes are continu- 
ally being updated and debugged. While many features of these codes 
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are similar, and so agreement between spectra tends to be good in generic 
parameter space regions, there are some differences as well. In particu- 
lar, the codes SuSpect , Softsusy and Spheno all adopt a sharp cut-off 
scale between the MSSM and SM effective theories. Allowance for the 
sharp cut-off is compensated for by log terms in the 1-loop sparticle and 
Higgs boson mass corrections. The ISASUGRA code instead adopts a "tower 
of effective theories" approach, and incorporates threshold corrections in 
the 1-loop RGEs. Here, the beta-functions changes each time a sparticle 
mass threshold is passed over. One loop corrections to non-mixing sparticle 
masses are implemented at each sparticle's mass scale, so all logs are min- 
imized. Sparticles that mix have all their SUSY parameters evaluated at 
the Msusy = ^/ m t L m t R scale due to a need for consistency amongst the 
various soft terms that enter the mass matrices. 17 In this way, better accu- 
racy is expected in cases where the sparticle mass spectra is spread across a 
large energy range, as happens- for instance- in focus point or split SUSY, 
where scalars are at multi-TeV values or beyond, whereas gauginos can be 
quite light. 



1.3. Sparticle production and decay codes 
1.3.1. Production cross sections 

The multi-purpose event generators Isajet ,Pythia ,Herwig , SUSYGeN 
and Sherpa all have a complete set of tree-level SUSY particle production 
reactions encoded, and can be used to calculate tree-level sparticle pro- 
duction cross sections. In the case of Pythia or Herwig , the LHA files 
from spectra generators can be used as input to calculate these, or general 
SUSY parameter inputs are allowed. Isajet does not allow LHA input 
since it has its own spectra generator. The Isajet code also calculates all 
sparticle and Higgs boson production reactions for e + e~ colliders including 
variable beam polarization, and bremsstrahlung and beamstrahhmg. 18 The 
Spheno code also calculates lowest order e + e~ — > SUSY cross sections. 

The code Prospino 19 has been created to calculate all 2 — > 2 super- 
symmetric production cross sections at hadron colliders at both leading 
order (LO) and next-to- leading order (NLO) in QCD. The current version 
of Prospino takes LHA files as its input format. 
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1.3.2. Decay widths and branching fractions 

The programs ISASUSY and ISASUGRA also calculate all sparticle and Higgs 
boson 1 — ► 2-body and 1 — > 3-body decay widths and branching fractions 
(BFs). These widths and BFs are output in Isajet standard output files, 
and are used internally for event generation. The chargino and neutralino 
branching fractions are sensitive to the parameter tan (3 in that at large 
tan/3, decays to third generation quarks and leptons are enhanced relative 
to decays to first/second generation fermions. 

The program SusyHIT 20 is a relatively new release that combines SuS- 
PECT with the branching fraction codes SDecay and HDecay to also gen- 
erate a table of sparticle and Higgs boson decay widths and branching 
fractions. Some of the decay modes in SusyHIT are calculated at NLO in 
QCD. 

The program Spheno also computes sparticle decay widths and branch- 
ing fractions. 

The branching fractions from Isajet , SusyHIT and Spheno all seem 
to enjoy excellent agreement with each other. The branching fractions of 
all these codes can be input to event generators via the LHA input /output 
files. Early versions of Herwig took branching fraction inputs from the 
IsaWIG output files. 

Care must be taken in extracting branching fractions for neutralinos 
and charginos computed internally from Pythia in that they may not be 
valid at large tan/3 values ~ 10, since Yukawa couplings, mixing effects, 
and decays through intermediate Higgs bosons are neglected. 

Some specialized codes are available for calculating decays modes of the 
SUSY Higgs bosons. These include FeynHiggs, 21 which calculates MSSM 
Higgs boson masses at two-loop level, along with branching fractions, CP- 
SUPERH 22 which calculates Higgs boson branching fractions including CP- 
violating parameters, and NHMDecay, 23 which calculates Higgs boson 
masses and branching fractions in the next-to-minimal supersymmetric 
Standard Model (NMSSM). 

1.4. Event generators 

Supersymmetric models can be used to calculate sparticle masses and mix- 
ings, which in turn allow for a prediction of various sparticle production 
rates and decay widths into final states containing quarks, leptons, pho- 
tons, gluons (and LSPs in imparity conserving models). However, quarks 
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and gluons are never directly measured in any collider detector. Instead, 
detectors measure tracks of (quasi-)stable charged particles and their mo- 
menta as they bend in a magnetic field. They also measure energy deposited 
in calorimeter cells by hadrons, charged leptons and photons. There is thus 
a gap between the predictions of supersymmetric models in terms of final 
states involving quarks, gluons, leptons and photons, with what is actually 
detected in the experimental apparatus. This gap is bridged by event gen- 
erator computer programs. Once a collider type and supersymmetric model 
is specified, the event generator program can produce a complete simulation 
of the sorts of scattering events that are to be expected. The final state of 
any simulated scattering event is composed of a listing of electrons, muons, 
photons and the long-lived hadrons (pions, kaons, nucleons etc.) and their 
associated 4-vectors that may be measured in a collider experiment. 

The underlying idea of SUSY event generator programs is that for a 
specified collider type (e + e~, pp, pp, • • • ) and center of mass energy, the 
event generator will, for any set of model parameters, generate various spar- 
ticle pair production events in the ratio of their production cross sections, 
and with distributions as given by their differential cross sections. More- 
over, the produced sparticles will undergo a (possibly multi-step cascade) 
decay 24 into a partonic final state, according to branching ratios as fixed 
by the model. Finally, this partonic final state is converted to one that is 
comprised of particles that are detected in an experimental apparatus. By 
generating a large number of "SUSY events" using these computer codes, 
the user can statistically simulate the various final states that are expected 
to be produced within the framework of any particular model. 

Several general purpose event generator programs that incorporate 
SUSY are currently available, including Isajet, 8 Pythia, 25 Herwig, 26 
SUSYGeN 27 and Sherpa. 28 These include usually just the leading order 
2^2 SUSY production processes. 

In addition, specific 2 — > n (n < 6) SUSY reactions may be generated by 
such programs as CompHEP , CalcHEP , MadGraph , SUSY-Grace 
, AmegicH — h and O'Mega . The output of these latter programs can be 
interfaced with the Pythia or Herwig programs to yield complete scat- 
tering event simulations by generating output in the Les Houches Event file 
(LHE) format. (Isajet generates LHE output, but does not accept LHE files 
as input, since it includes its own mass and branching fraction generator). 
Ideally, event generator programs should be flexible enough to enable simu- 
lation of SUSY events from a variety of models such as mSUGRA, GMSB, 
AMSB etc.. This is usually accomplished nowadays by reading in the LHA 
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model files into the event generators Pythia or Herwig . Since Isajet 
does its own spectra calculation, it only outputs LHA files, but does not 
accept them as input. 

The simulation of hadron collider scattering events may be broken up 
into several steps, as illustrated in Fig. 11.11 The steps include: 

• the perturbative calculation of the hard scattering subprocess in the 
parton model, and convolution with parton distribution functions 
(PDFs), 

• inclusion of sparticle cascade decays, 

• implementation of perturbative parton showers for initial and final 
state colored particles, and for other colored particles which may 
be produced as decay products of heavier objects, 

• implementation of a hadronization model which describes the for- 
mation of mesons and baryons from quarks and gluons. Also, un- 
stable particles must be decayed to the (quasi-)stable daughters 
that are ultimately detected in the apparatus, with rates and dis- 
tributions in accord with their measured or predicted values. 

• Finally, the debris from the colored remnants of the initial beams 
must be modeled to obtain a valid description of physics in the 
forward regions of the collider detector. 

For e + e~ collider simulations, in addition we have to allow for the possibility 
of polarized initial beams, and beam-strahlung effects. 

1.4.1. Hard scattering 

The hard scattering and convolution with PDFs forms the central calcula- 
tion of event generator programs. The calculations are usually performed 
at lowest order in perturbation theory, so that the hard scattering is either 
a2— >2or2^1 scattering process. Matrix element generators are usually 
used for 2 — > n processes, with n > 3. 

For supersymmetric particle production at a high energy hadron collider 
such as the LHC, a large number of hard scattering subprocesses are likely to 
be kinematically accessible. Each subprocess reaction must be convoluted 
with PDFs so that a total hadronic cross section for each reaction may be 
determined. The <5 2 -dependent PDFs commonly used are constructed to 
be solutions of the Dokshitzer, Gribov, Lipatov, Altarelli, Parisi (DGLAP) 
QCD evolution equations, which account for multiple collinear emissions of 
quarks and gluons from the initial state in the leading log approximation. 
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Event generation in LL - QCD 

1) Hard scattering / convolution with PDFs 

2) Intial / final state showers 

3) Cascade decays 

4) Hadronization 

5) Beam remnants 

\ / 

Fig. 1.1. Steps in any event generation procedure. 

As Q 2 increases, more gluons are radiated, so that the distributions soften 
for large values of x, and correspondingly increase at small x values. Use 
of a running QCD coupling constant makes the entire calculation valid at 
leading log level. 

Once the total cross sections are evaluated for all the allowed subpro- 
cesses, then reactions may be selected probabilistically (with an assigned 
weight) using a random number generator. This will yield sparticle events 
in the ratio predicted by the particular model being simulated. 

For sparticle production at e + e~ colliders, it may also be necessary to 
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convolute with electron and positron PDFs to incorporate bremsstrahlung 
and beamstrahlung effects. In addition, if beam polarization is used, then 
each subprocess cross section will depend on beam polarization parameters 
as well. 

1.4.2. Parton showers 

For reactions occurring at both hadron and lepton colliders, to obtain a 
realistic portrait of supersymmetric (or Standard Model) events, it is nec- 
essary to account for multiple non-collinear QCD radiation effects. The 
evaluation of the cross section using matrix elements for multi-parton fi- 
nal states is prohibitively difficult. Instead, these multiple emissions are 
approximately included in an event simulation via a parton shower (PS) 
algorithm. They give rise to effects such as jet broadening, radiation in the 
forward regions and energy flow into detector regions that are not described 
by calculations with only a limited number of final state partons. 

In leading log approximation (LLA), the cross section for single gluon 
emission from a quark line is given by 

da = a ^-^P qq (z)dz, (1.1) 

where <7o is the overall hard scattering cross section, t is the intermedi- 
ate state virtual quark mass, and P qq (z) — | ^ X ^-z ) comc ides with the 
Altarelli-Parisi splitting function for q' — » qg for the fractional momentum 
of the final quark z = j^jj < 1. Interference between various multiple 
gluon emission Feynman graphs, where the gluons are ordered differently, 
is a subleading effect. Thus, Eq. (|l.l|l can be applied successively, and gives 
a factorized probability for each gluon emission. The idea behind the PS 
algorithm is then to use these approximate emission probabilities (which 
are exact in the collinear limit), along with exact (non-collinear) kinemat- 
ics to construct a program which describes multiple non-collinear parton 
emissions. Notice, however, that the cross section (|1.1|) is singular as t — > 
and as z — ► 1, i.e. in the regime of collinear and also soft gluon emission. 
These singularities can be regulated by introducing physically appropriate 
cut-offs. A cutoff on the value of \t\ of order \t c \ ~ 1 GeV corresponds to 
the scale below which QCD perturbation theory is no longer valid. A cutoff 
on z is also necessary, and physically corresponds to the limit beyond which 
the gluon is too soft to be resolved. 

The PS algorithms available vary in their degree of sophistication. The 
simplest algorithm was created by Fox and Wolfram in 1979. 29 Their 
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method was improved to account for interference effects in the angle-ordered 
algorithm of Marchesini and Webber. 30 In addition, parton emission from 
heavy particles results in a dead-cone effect, where emissions in the direc- 
tion of the heavy particle are suppressed. Furthermore, it is possible to 
include spin correlations in the PS algorithm. 

PS algorithms are also applied to the initial state partons. In this 
backwards shower algorithm is most efficient, which develops the 
emissions from the hard scattering backwards in time towards the initial 
state. The backward shower algorithm developed by Sjostrand 31 makes use 
of the PDFs evaluated at different energy scales to calculate the initial state 
parton emission probabilities. 

1.4.3. Cascade decays 

Not only are there many reactions available via which SUSY particles may 
be produced at colliders, but once produced, there exist many ways in which 
superparticles may decay. For the next-to- lightest SUSY particle (NLSP), 
there may be only one or at most a few ways to decay to the LSP. Thus, 
for a collider such as LEP or even the Fcrmilab Tevatron, where only the 
lightest sparticles will have significant production rates, we might expect 
that their associated decay patterns will be relatively simple. However, the 
number of possible final states increases rapidly if squarks and gluinos that 
can decay into the heavier charginos and neutralinos are accessible, and 
the book-keeping becomes correspondingly more complicated. Indeed, at 
the CERN LHC, where the massive strongly interacting sparticles such as 
gluinos and squarks are expected to be produced at large rates, sparticle 
cascade decay patterns can be very complex. As a rough estimate, of order 
10 3 subprocess cross sections may be active at LHC energies, with of order 
10 decay modes for each sparticle. Naively, this would give of order 10 5 
2 — > n-body subprocesses that would need to be calculated. 

Monte Carlo event generators immensely facilitate the analysis of signals 
from such complex cascade decays, especially in the case where no single 
decay chain dominates. An event generator can select different cascade 
decay branches by generating a random number which picks out a particular 
decay choice, with a weight proportional to the corresponding branching 
fraction, at each step of the cascade decay. Quarks and gluons produced 
as the end products of cascade decays will shower off still more quarks and 
gluons, with probabilities determined by the PS algorithm. 

The procedure that we have just described is exact for cascade decays 
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of spinless particles into two other spinless particles at each step in the 
cascade. This is because the squared matrix element is just a constant, 
and there are no spin correlations possible. This is not true in general 
and in many cases, it can be very important to include the decay matrix 
element and/or spin correlations in the calculation of cascade decays of 
sparticles. A general method for incorporating spin correlations based on 
density matrices has been put forth by Richardson, 32 and incorporated into 
Herwig . 

Spin correlation effects are especially important for precision measure- 
ments at e+e~ linear colliders. While retaining spin correlations may be 
less crucial in many situations at a hadron collider, this is not always the 
case. For instance, relativistic t~ leptons produced from W decay are al- 
ways left-handed, while those produced from a charged Higgs decay are 
always right-handed. Likewise, the polarization of the taus from f\ decays 
depends on the stau mixing angle. Since the undetectable energy carried 
off by v T from tau decay depends sensitively on the parent tau hclicity, 
it is necessary to include effects of tau polarization in any consideration 
involving the energy of "tau jets" . By evaluating the mean polarization of 
taus in any particular process, these effects can be incorporated, at least 
on average, into event generator programs such as Isajet . Of course, such 
a procedure would not include correlations between decay products of two 
taus produced in the same reaction. 

Another aspect is to include appropriately the complete 3-body decay 
matrix elements. While some programs merely use a flat phase space distri- 
bution, Isajet and Herwig include pre-programmed exact decay matrix 
elements. 

1.4.4. Models of hadronization 

Once sparticles have been produced and have decayed through their cas- 
cades, and parton showers have been evolved up to the point where the par- 
tons have virtuality smaller than ~ 1 GeV 2 , the partons must be converted 
to hadrons. This is a non-perturbative process, and one must appeal to phc- 
nomenological models for its description. The independent hadronization 
(IH) model of Field and Feynman 33 is the simplest such model to imple- 
ment. In this picture, a new quark anti-quark pair qiqi can be created in 
the color field of the parent quark go- Then the go9i pair can turn into a 
meson with a longitudinal momentum fraction described by a phenomeno- 
logical function, with the remainder of the longitudinal momentum carried 
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by the quark q x . This process is repeated by the creation of a <72<?2 pair in 
the color field of qi, and so on down the line to q n q n . A host of mesons 
are thus produced, and decayed to the quasi-stable tt, K, ■ ■ ■ mesons ac- 
cording to their experimental properties. The final residual quark q n will 
have very little energy, and can be discarded without significantly affecting 
jet physics. Finally, a small transverse momentum can be added according 
to a pre-assigned Gaussian probability distribution to obtain a better de- 
scription of the data. Quark fragmentation into baryons is also possible by 
creation of diquark pairs in its color field, and can be incorporated. The 
IH scheme, with many parameters tuned to fit the data, will thus describe 
the bulk features of hadronization needed for event simulation programs. 

The string model of hadronization developed at Lund 34 is a more so- 
phisticated model than IH, which treats hadron production as a universal 
process independent of the environment of the fragmenting quark. In the 
string model, a produced quark-antiquark pair is assumed to be connected 
by a color flux tube or string. As the quark-antiquark pair moves apart, 
more and more energy is stored in the string until it is energetically favor- 
able for the string to break, creating a new quark-antiquark pair. Gluons 
are regarded as kinks in the string. The string model correctly accounts 
for color flow in the hadronization process, as opposed to the IH model. 
In e + e~ — ► qqg (3-jet) events, the string model predicts fewer produced 
hadrons in the regions between jets than the IH model, in accord with 
observation. 

A third scheme for hadronization is known as the cluster hadronization 
model. 35 In this case, color flow is still accounted for, but quarks and anti- 
quarks that are nearby in phase space will form a cluster, and will hadronize 
according to preassigned probabilities. This model avoids non-locality prob- 
lems associated with the string hadronization model, where quarks and 
antiquarks separated by spacelike distances can affect the hadronization 
process. 

1.4.5. Beam remnants 

Finally, at a hadron collider the colored remnants of the nucleon that did 
not participate in the hard scattering must be accounted for. These beam 
remnant effects produce additional energy flow, especially in the far forward 
regions of the detector. A variety of approaches are available to describe 
these non-perturbative processes, including models involving Pomeron ex- 
change and multiple scatterings. In addition, the beam remnants must 
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be hadronized as well, and appear to require a different parametrization 
from "minimum bias" events where there are only beam jets but no hard 
scattering. 

1.4.6. Multi-purpose event generators 

Publicly available event generators for SUSY processes include, 

• ISAJET: (H. Baer, F. Paige, S. Protopopescu and X. Tata), 
http : //www . nhn . ou . edu/~isaj et/ 

• PYTHIA: (T. Sjostrand, L. Lonnblad and S. Mrenna), 
http : / /www . thep . lu . se/~torb j orn/Pythia . html 

• HERWIG: (G. Corcella, I. G. Knowles, G. Marchesini, S. Moretti, 
K. Odagiri, P. Richardson, M. Seymour and B. R. Webber), 
http : / /hepwww . rl . ac . uk/theory/ seymour /herwig/ 

• SUSYGEN: (N. Ghodbane, S. Katsanevas, P. Morawitz and E. 
Perez), 

http : / /lyoinf o . in2p3 . f r/ susygen/susygen3 .html 

• SHERPA: (T. Gleisberg, S. H6chc, F. 
Krauss, M. Schonherr, S. Schumann, F. Siegert and j. Winter) 
http : //projects .hepf orge . org/ sherpa/dokuwiki/doku .php 

The event generator program Isajet was originally developed in the late 
1970's to describe scattering events at the ill-fated ISABELLE pp collider 
at Brookhaven National Laboratory. It was developed by F. Paige and 
S. Protopopescu to generate SM and beyond scattering events at hadron 
and e + e~ colliders. H. Baer and X. Tata, in collaboration with Paige 
and Protopopescu, developed Isajet to give a realistic portrayal of SUSY 
scattering events. Isajet uses the IH model for hadronization, and the 
original Fox- Wolfram (Sjostrand) PS shower algorithm for final state (initial 
state) parton showers. It includes an n-cut Pomeron model to describe 
beam-jet evolution. 

The event generator Pythia was developed mainly by T. Sjostrand in 
the early 1980s to implement the Lund string model for event generation. 
Pythia uses the FW virtuality-ordcred shower model, but with an angle- 
ordered veto. S. Mrenna contributed the inclusion of SUSY processes in 
Pythia . 

The event generator Herwig was developed in the mid-1980s to describe 
scattering events with angle-ordered parton showers, which accounted for 
interference effects neglected in the FW shower approach. Herwig implc- 
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ments a cluster hadronization model. Herwig is notable in that it includes 
sparticle production and decay spin correlations using density matrix tech- 
mques. 

The program SUSYGeN was developed by S. Katsanevas and P. 
Morawitz to generate e + e~ — > SUSY events for the LEP experiments. 
SUSYGeN interfaces with Pythia for hadronization and showering. 
SUSYGeN has since been upgraded to also generate events for hadron 
colliders. 

The program Sherpa was developed as a new generation event gen- 
erator in the C + + language. It calculates subprocess reactions using 
Amegic++ . It includes its own shower and cluster hadronization rou- 
tines. 

1.4.7. Matrix element generators 

For generating various 2 — > n scattering reactions using complete matrix 
elements, a number of automated tree-level codes are available. 

The code CompHEP by E. Boos et al. 36 is designed to take one di- 
rectly from a Lagrangian to distributions. Feynman rules can be calculated 
using the LanHEP code, 37 and then CompHEP will generate the squared 
matrix element by constructing the squared amplitude, taking traces, and 
storing the output as subroutines. CompHEP also includes code for doing 
the phase space integration, convolution with PDFs, and after integration, 
numerical output, or output in terms of histograms. 

The code CalcHEP 38 by Pukhov, Belyaev and Christensen is very 
similar to CompHEP , and was in fact created as a spin-off by some of the 
original authors of CompHEP . 

The code MadGraph /MadEvent was developed by Stelzer and Long 
and others. 39 It allows the user to input initial and final state particles, 
and then generates all Feynman diagrams along with a subroutine which 
evaluates the scattering amplitude as a complex number using the Helas 
hclicity amplitude subroutines developed by Hagiwara and Murayama. 40 
Since MadGraph directly evaluates the amplitude, and not amplitude 
squared, computational sampling of the squared matrix element should be 
faster than programs which evaluate traces over gamma matrices. The 
latest versions of MadGraph , updated to MadEvent , will convolute 
with PDFs and perform phase space integration and evaluate distributions 
as well. A number of models for BSM physics, including the MSSM , are 
available in MadGraph /MadEvent . 
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The program O'Mega by Ohl, Reuter and Schwinn, also generates 
tree-level SM and MSSM amplitudes, and can work in concert with the 
Whizard program for event generation. 41,42 

The program SUSY-Grace by Tanaka, Kuroda, Kaneko, Jimbo and 
Kon also generates SM and MSSM amplitudes, and generates scattering 
events in association with the Grappa program. 43,44 

1.4.8. Les Houches Event (LHE) files 

A Les Houches Event (LHE) file format has been proposed 45 which allows 
for a simple communication between parton level event generators and all 
purpose generators such as Pythia and Herwig . This is particularly 
useful when matrix element generators like CalcHEP or MadGraph are 
used, but the user needs a complete event output including parton showers, 
hadronization and underlying event simulation. 

The LHE file is an ascii file which includes lines pertaining to the 
generator initialization. It then follows with a listing of partons (particle ID 
code), their associated 4- vectors and color flow information. The generators 
Pythia and Herwig then can read in these files, to add on showering, 
hadronization and underlying event. A sample SUSY event in LHE format 
is listed below. It lists a reaction with sg — -> sg, with s — > sxl and g — > 
dd and then d — > d\\. After listing a line of event characteristics, the 
event listing follows. The first column corresponds to particle ID, 2nd 
column to stability of particle, 3rd and 4th columns list the source of the 
particle, 5th and 6th columns relate to color flow, and 7th column is the 
re-component of the energy-momentum four-vector. The four-vector listing 
has been truncated to fit on the page. 

<event> 

10 2160 

3 -1 

21 -1 

2000003 2 

1000021 2 

1000022 1 
3 1 

2000002 2 

-2 1 

1000022 1 

2 1 

</event> 
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102 



.768145E+06 

0.000000E+00 

101 0.000000E+00 

0.220402E+03 

103 -.220402E+03 

0.778961E+02 

0.142506E+03 

-.185972E+03 

103 -.344299E+02 

. 800434E+02 

-.266016E+03 
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1.5. Dark matter codes 

In response to the increasing precision of data corresponding to the density 
of dark matter in the universe, several public codes have been developed 
which evaluate key astrophysical observables in supersymmetric (and other) 
models. 



1.5.1. DarkSUSY 

The DarkSUSY code, developed by Gondolo et al., 46 evaluates the relic 
density of neutralino dark matter in SUSY models. DarkSUSY computes all 
relevant neutralino annihilation and co-annihilation processes, and solves 
the Boltzmann equation to output the current density of neutralino CDM. It 
accepts input files from Isajet/Isasugra or from LHA input files. DarkSUSY 
also calculates: spin-independent and spin-dependent ncutralino-nucleon 
scattering rates (direct WIMP detection) , and indirect neutralino detection 
rates, such as: muon flux from neutralino annihilation in the core of earth 
or sun, flux of 7 rays, pa, e + s and ds from neutralino annihilation in the 
galactic core or halo. The halo annihilation rates all depend on an assumed 
form for the galactic dark matter density profile. 

1.5.2. Micromegas 

MicrOMEGAS was developed by Belanger et al., 47 and also evaluates the 
neutralino relic density due to all annihilation and co-annihilation processes. 
It also computes the WIMP relic density for a variety of other non-SUSY 
models. It also outputs neutralino direct and indirect detection rates, b — > 
S7 branching fraction and (g — 2) p . 

1.5.3. Isatools 

ISATOOLS is part of the Isajet package. It includes a subroutine IsaReD 48 
to evaluate the neutralino relic density, the direct neutralino detection rates 
via spin-independent and spin-dependent scattering, the b — > s-f branching 
fraction, (g — 2)^ USY , BF(B S — > an d the thermally averaged neu- 

tralino annihilation cross section, which is key input to neutralino halo 
annihilation calculations. 
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1.6. Parameter fitting codes 

If supersymmetry is indeed discovered at the Tevatron, LHC and/or a linear 
e + e~ collider, then an exciting task will be to make precision measurements 
of all sparticle masses, spins, couplings and mixings. Once these are known, 
then, if the MSSM is indeed the correct effective theory all the way from 
M wea k to Mqut, it is possible to map out the GUT scale values of the soft 
SUSY breaking parameters. Once these are known, important information 
will be gained which will allow for the construction of SUSY models at or 
beyond the GUT scale. Two such codes are available which accomplish this 
task: Sfitter 49 and Fittino. 50 

1.7. SPA convention 

The supersymmetry parameter analysis (SPA) project 51 is an attempt to 
achieve co-ordination between the various sparticle mass generation codes, 
event generators, relic density codes, and parameter fitting codes, with a 
goal in mind to determine the fundamental SUSY Lagrangian. In the SPA 
convention, all programs should input /ouput SUSY parameters in the DR 
scheme at the Q — 1 TeV scale. Once this benchmark is set, then all 
remaining calculations may proceed from this common agreed upon point. 

1.8. Summary 

In the past decade, there has been an explosion of interest in supersymmetry 
phenomenology. This is exhibited in part by the corresponding development 
of numerous computational tools to aid in supersymmetry calculations for 
expected collider events and for dark matter observables. Supersymmetry 
has certainly been an enduring theme in high energy physics. Hopefully, 
at the dawn of the LHC era, we are on the verge of actual discovery of 
supersymmetry. In this case, many of these tools for SUSY will be put to 
good hard use, as the community analyzes the upcoming collider data. 

We expect that new tools for SUSY will emerge, which will be more fo- 
cused on the new matter states that might appear. As an example, if SUSY 
is discovered, then the MSSM (or perhaps NMSSM) may become the new 
SM, and radiative corrections will have to be calculated for any remaining 
production and decay reactions, and in a form suitable for embedding in 
event generator programs. The clues we find pertaining to dark matter 
will impact on all astrophysical codes. In addition, new tools should also 
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emerge that facilitate model building, as the clues we expect to emerge from 
the data point the way to a new paradigm in physics beyond the Standard 
Model. 
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